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Abstract 

The method of fundamental solution (MFS) is an efficient meshless method 
for solving a boundary value problem in an exterior unbounded domain. The 
numerical solution obtained by the MFS is accurate, while the corresponding 
matrix equation is ill-conditioned. A modified MFS (MMFS) with the proper 
basis functions is proposed by the introduction of the modified Trcfftz method 
(MTM). The concrete expressions of the corresponding condition numbers and 
the solvability by these methods are mathematically proven. Thereby, the opti- 
mal parameter minimizing the condition number is also mathematically given. 
Numerical experiments show that the condition numbers of the matrices corre- 
sponding to the MTM and the MMFS are reduced and that the numerical solu- 
tion by the MMFS is more accurate than the one by the conventional method. 

Keywords: Exterior unbounded domain, Condition number, Laplace 
equation, Method of fundamental solution, Modified Trefftz method, 



1. Introduction 

The method of fundamental solutions (MFS) is a truly meshless numerical 
method for easily and rapidly solving boundary value problems of elliptic type. 
The approximate solution by the MFS is expressed by a linear combination of 
fundamental solutions to a partial differential equation. Hence, the approximate 
solution automatically satisfies the partial differential equation, and it remains 
to consider boundary conditions. In this sense, the MFS is a boundary method 
like the boundary element method (BEM). The BEM requires singular inte- 
grals, while the MFS does not require any treatments for the singularity of the 
fundamental solutions, which is an advantage of the MFS. 
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Mathon and Johnston [8] first obtained numerical solutions by applying the 
MFS. Bogomolny [1] proved the convergence property. Katsurada and Okamoto 
[4], [6] proved uniqueness of numerical solutions and the convergence property 
and mathematically discussed the suitable distribution of source points. 

Another advantage of the MFS is to directly solve exterior problems, which 
is different from the finite element method (FEM) and other meshless methods 
using radial basis functions. Actually, the MFS can be applied directly to the 
exterior problems to obtain accurate solutions. Even though the basis functions 
used in the MFS do not satisfy a condition at infinity, the accuracy is not bad in 
the whole computational domain. Katsurada [5] proposed that the approximate 
solution for the exterior problems should be defined by a linear combination 
of the proper basis functions satisfying the governing equation as well as the 
condition at infinity. Although many researchers seem to still use the MFS with 
the conventional basis functions, we should use the MFS with the proper basis 
functions to obtain more accurate solutions in the whole domain. 

On the other hand, the Trefftz method is also known as a numerical meshless 
method for solving boundary value problems. The approximate solution by the 
Trefftz method is expressed by a linear combination of functions satisfying the 
governing equation. Hence, similar to the MFS, it is sufficient if the boundary 
condition is considered. This method is also regarded as a boundary method. 

It is well known that the coefficient matrices derived from the MFS and 
the Trefftz method are highly ill-conditioned. Numerical solutions to the ill- 
conditioned matrix equations are unstable. We need to improve the ill-conditioning. 

Ramachandran [9] applied the singular value decomposition (SVD) to the 
MFS to improve the accuracy. Chen et al. [2], [3] discussed the equivalence 
between the Trefftz method and the MFS. Liu [7] proposed the modified MFS 
(MMFS) for interior problems by introducing the modified Trefftz method (MTM). 
Both of the MTM and the MMFS can drastically reduce the condition numbers 
of the corresponding matrices. 

In this paper, we propose a modified MFS for solving the exterior problem 
according to the papers cited above. The solution of the exterior problem defined 
in Section 2 is discretized by the MFS with the proper basis functions satisfying 
the condition at infinity in Section 3. Section 4 proposes a modified MFS based 
on the modified Trefftz method with a characteristic length, which plays an 
important role for the condition number. Section 5 shows some mathematical 
results about the condition numbers and the uniqueness of the approximate 
solutions by the MTM and MMFS. Thereby, the optimal characteristic length 
that minimizes the condition number corresponding to the MTM can explicitly 
be given in a mathematical form. Section 6 shows the effectiveness of the MTM 
and the MMFS through numerical experiments. 

2. Problem Setting 

Let fio be a two dimensional bounded domain enclosed by the boundary 
r = {(r, 9) : r = p(9), < 9 < 2ir} in the polar representation, where p(9) is 
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a smooth function. Let Q := R 2 \ fig be the unbounded domain outside the 
domain fio- Then, we consider the following exterior Dirichlct problem of the 
Laplace equation: find u e C 2 (0) n C(fi) such that 

Au = in Ct, (1) 

u = f on L, (2) 

u(x) = 0(|a:| _1 ) as |x| — »■ oo, (3) 

where / is a continuous function defined on L. 

From now on, we identify the Cartesian plane R 2 with the complex plane 
C. We denote the (j, k) component of a matrix Q by Qj t k- 

3. Discretization by the Method of Fundamental Solutions (MFS) 

3.1. The conventional MFS with the conventional basis functions 

The fundamental solution of the Laplace equation in two dimensions is de- 
fined as 

G*(r) :=-^lnr 

for r = \z\ = yj x 2 + y 2 , which is a solution to 

—AG*(r) = S(r) 

with the Dirac delta distribution S. 

We distribute the source points {Cj}|Li along a circle outside the domain fl. 
The basis functions are defined as 

G j (z):=-2nG*(\z-( j \)=ln\z-C j \; {C,}f=i C Tf = fi . (4) 

Then, the exact solution u can be approximated by a linear combination of the 
basis functions as follows: 

JV 

u(z) !=a un(z) := y] WjGj(z), Vz€fi, (5) 
i=i 

where C i? are expansion coefficients to be determined below. Since 

the basis functions (4) have no singular points in 0, the approximate solution 
un satisfies the Laplace equation (1) in the domain f2. It remains to consider the 
boundary condition to find {wj}^^. We use the boundary collocation method 
since it is impossible that u^ exactly satisfies the boundary condition. We 
distribute the collocation points {zk}^ =1 on the boundary T. Substituting (5) 
at the collocation points into (2), we have 

JV 

^wjGjizk) = f(z k ), fc = i,2,...,iV; Wf =1 cr 

3=1 
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or in the matrix form: 

Aw = /, (6) 

where the matrix A = (Akj) e R NxN and the vectors w = (wj) e R N , 
f = (f k ) € R N are defined by 

A kJ -Gjizk), fc = l,2,...,JV; j = 1, 2, . . . , N, 
A :=/(**)> fc = l,2,...,7V. 

If ^4 is not singular, we can solve (6) and obtain the approximate solution un 
by substituting w into (5). 

For any point z = re 10 e 17 and the source point Q — Re l ^ j € f2o with the 
imaginary unit i — \J — 1, we have 

k-Cj| 2 = K 9 -i?e^| 2 = r 2 +i? 2 -2ri?cos(<9-<^) 
= 0(r 2 ), r — > oo, 

from which we know 

Gj(z) = In \z - Cj| = O(lnr), r -> oo. 

Hence, we obtain 

AT 

ujv(z) = ^ WjGj(z) = O(lnr), r -> oo, (7) 

which does not satisfy the condition (3). Therefore, the basis functions (4) are 
not proper for the exterior problem. 

3.2. The conventional MFS with modified basis functions 

We need to use another basis functions to approximate a solution that con- 
verges to zero at infinity. We define the following modified basis functions [5] : 

Gj{z) = In \ Z ~ ^ = In \z - Q\ - In \z\, (8) 
\z\ 

where Gj(z) satisfies AGj = in O since both of ln|z — Q\ and ln|z| are 
the fundamental solutions. For any point z = re 10 e il and the source point 
(j = Re l ^> e O , we have 

|z-C 3 | 2 _ r 2 + R 2 -2rRcos{6~(j) j ) 
\z\ 2 r 2 

= 1+ (*) 2 " 2 (f) «»(*-^) 
= l + 0(r" 1 ), r^oo, 
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from which we know 

Gj(z) = In ^ 7,^ = In y/l + Oir- 1 ) = O^ 1 ), r -> oo. (9) 
Hence, we obtain 

AT 

ujv(z) = ^ WjGj(z) = 0(r _1 ), r — > oo, 
i=i 

which satisfies the condition (3). Therefore, we can see that (8) are the proper 
basis functions for solving the exterior problem. In this paper, the approximate 
solution 

N 

u N (z) = ^2w j G j (z), VzeO (10) 
i=i 

is called the MFS with the modified basis functions. When we use this method, 
in (6) we replace the matrix A into A defined by A k j := Gj(zk)- 

4. A modified MFS with modified basis functions 

4--1- The modified Trefftz method (MTM) 

The approximate solution by the modified Trefftz method (MTM) [7] is given 

by 

u M (r,9) := a a + ^2 (a k cos k6 + b k sin k8) I — ) , r>p(6), < 6 < 2tt, 

k=l \ r / 

(11) 

where the characteristic length R is often taken as 

Ro < Pmin := min p(6). (12) 

This method is reduced to the conventional Trefftz method if Rq = 1. We can 
see that (3) implies that an = 0. But, we treat ao as an unknown coefficient for 
reasons of expediency (see Remark 1). Using the boundary collocation method, 
we have 

m ( R \ k 

u M (p 3 ,0j) = a + ^2(a k cosk6j +b k smke 3 ) ( — ) = g h j = 1, 2, . . . , N, 
k=i \ p i ' 

(13) 

where 

Pj := ptfi), 9j := f(6j), Oj := 27r(j - l)/N. 
Therefore, we obtain the matrix equation: 

Sy = g, (14) 
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where 



y := (a ,ai,6i, . . . ,a M ,b M ) T € R 2M+1 , 
and S € ^ X(2M+1) is denned by 



9 := (5i,52,---,3at) T € 



AT 



ft 



coskOj, Sj^k+i ■= 



Pi 



sin 



(j = 1,2,. ..,7V; fc = l,2,...,M). 



When Ro satisfies (12), every element of S is less than or equal to 1. We take 
N = 2M + 1 to make S a square matrix. 

Let Zj — pje te i in (6). Then, / coincides with g. 

Remark 1. Since we know that ao = 0, the number of unknown coefficients 
are 2M in essentials. We can remove the first column of S, which we denote by 



S' . If we consider the square matrix S' € R 
see that 



2MxN 



with N — 2M, then we can 



2nd — 1) 

sin MOj = sinM — ^ — - = smn{j - 1) = (j = 1,2,..., N), 

which implies Sj 2 m ~ f° r 3 ~ 1, 2, . . . , N. Hence, S' becomes singular. 
Although we can avoid the singularity of S' by taking 9j = 2ir(j — 1)/(N + 1), 
we treat ao as an unknown coefficient to use S in this paper. 



4.1.1. Transformation of the MFS to the MTM 

The modified basis function for the source point £ = Re 1 ^ can be transformed 
as follows: 



In- 



CI 



In 



re' 



Re^l 



In 



re* 



1 _ ^e^-S) 



Rohi ( 1 - -e* 

r 



1 /i? 



= -Re£ 



1 /i? 



fe=l 



fe=l 

oo 



o ife(0-e) = _ 



= -E 



i /i? 



cosfc(0 — (9) 



00 1 / R\ k 

fc=l V 7 



Hence, we can write the approximate solution by the MFS in the form: 

AT 



JV 

E 



° 1 / R • ; 



E 

fc=i 



cosk{6-6j) 



(15) 
(16) 



G 



where we put Q = Re 10 ' with Oj = 2w(j - l)/N (j = 1,2,..., N). If we use 
the conventional basis functions, the approximate solution can be written in the 
form: 



N 



Un{t, 6) — Wj In \z — (. 

N 



1/7? 



fe=l 



(17) 
(18) 



We know from (7) that (17) or (18) does not satisfy (3). But, if we impose the 
condition 

N 

^ Wj In r — (r^l) 



or 



E W J = °' 



(19) 



then (18) coincides with (16). Hence, we know that (18) satisfies (3) as long as 
(19) is imposed. 

On the boundary 

T = {(r,0):r = p{0), < 9 < 2ir}, 
the approximate solution (16) is written as 



Wj 



1 / R \ 



, s 1 ( R N k 



E E /,• V/-'; 



3 

oo AT 



cosk{6-6j) 
(cos cos + sin k9 sin fcfy) 



= -EE^- 



1 / R 



k=ij=i 



k \R 



(cos k9j cos fc# + sin k9j sin 



I? 



iV 



Truncating the infinite series X^feLi m (20) into a finite series X^fcLr 

m n i / R\ k 
un(p(9), 9) w — ^ ^ w 3 — I — I (cos cos + sin fc^- sin k0) 



(20) 



(21) 
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and comparing (21) with (11), we obtain 
1 / R \ k N 

a k = -- f —J cos kOj, fc = l,2,...,M, (22) 

3=1 
k N 

6^ = —j- ( ) y^WjSinfc^-, fc = 1,2,...,M. (23) 



1 / R\ 

(^J ^WjSinfc^-, fc = l,2,...,M. 



Adding the condition (19) to (22) and (23), we obtain the following matrix 
equation: 

y = Kw, (24) 
where K E i?( 2 ^+i)xiv ig dcfincd by 

1 / R \ X / R \ 

K i,j '■= !' ^2fej := -- ( --J cosfcflj, K 2k+i,j ■= -j ( -- J sinfctf, 

(fc = l,2,...,M; i = 1,2,. ..,7V). 

Substituting (24) into (14), we obtain the matrix equation 

SKw = f. (25) 

Remark 2. When we consider (16), it is not necessary to impose (19) in es- 
sentials. In addition, we know that ao = 0. Hence, we can remove the first row 
of K, which is denoted by K' e R 2MxN . Then, the square matrix K' with 
N = 2M becomes singular because of the same reason described in Remark 1. 
Therefore, we impose the condition (19) to use K. 

The matrix K for i? = 1, denoted by K\, becomes 

K ltj = 1, K 2 h,j = ~T Rk coskOj, K 2k +i.j = -~r Rk sinkOj 

(fc = l,2,...,M; j = l,2,...,N), 

while K for i? = R, denoted by K 2 , becomes 



-j-cosfctfj, K 2k +i,j = 

(fc = l,2,...,M; j = l,2,...,A0, 



Kij = 1, K 2k ,j = ~r cos fc6»j, K 2k+ ij = -- sink6j 



which is independent of R. 

After obtaining the unknown vector w, we substitute w into (15) not (21). 
We call this method the modified MFS (MMFS). Even if we substitute w into 
(21), the result should be the same as the result when we use (13). 

Remark 3. We can see that SK <E R NxN is always a square matrix for any M. 
Since we would like to find unknown w from the known /, it is not necessary 
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that S and K arc square matrices. However, we see from numerical experiments 
that a better approximate solution can be obtained when the matrices S and 
K are square. This numerical result is not shown in this paper. 

We notice that the matrix SK converges to the matrix A corresponding to 
the conventional MFS as M — > oo. In this sense, SK can be regarded as an 
approximation to A. 



5. Some Mathematical Results 

In this section, we show some theoretical results about the matrices S and 
K defined in the previous section. In particular, we obtain the optimal i? 
in a mathematical form. We prove the uniqueness of the numerical solutions 
obtained by the MTM and the MMFS. 

The following Propositions 1 and 2 have already been shown in [2] and [7]. 
In this paper, we give the complete and more simple proof of Proposition 2. 

Proposition 1. The matrix K can be decomposed into 

K = T R T e , 

where T R e fi (2M+1)x(2M+1) and T e E R^ M +^ N are defined as 

(T R ) hl = l, {T R ) kJ =0 (kjtj) 

1 / R N ' 



(T R ) 2k ,2k = (T fl )2fc+i,2fc+i = - ^ ^ ) (k = 1, 2, . . . , M) 



{Te)i,i = l (j = l,2,...,N), 
{T ) 2 k,j =cosk6 j , (Te) 2 k+i,j =smk9 J (k = 1, 2, . . . , M; j = 1, 2, . . . , AO- 
Proposition 2. The matrix Tg has the orthogonal property: 

TgT T = E. diag(2; 1, . . . , 1) g R (2M+l)x(2M+l)_ (2g) 

When N = 2M + 1, the inverse matrix Tg 1 can be written as 
(Z*- 1 )*! = (j = l,2,...,N), 

{Tg 1 ) j , 2k = lcoakO j , (T ~ 1 )j£k+i = sin kOj (k = 1, 2, . . . , M; j = 1, 2, . . . , N). 
Proof. Let w := cxp(27ri/iV). Then, we see exp(ikdi) = cv k< - l ^ 1 \ We notice that 



f, )S (l-i) = f JV (s = miV; m = 0,±l,±2,...) 
^ " \ {s^mN; m = 0, ±1, ±2, . . .) 



a; 
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since for s ^ mN, we can divide the following equation by u s — 1 ^ 0: 
(lo s - 1) ({oj 3 )"- 1 + (lo s ) n - 2 + ■ ■ ■ + lj s + 1) = (lo s ) n -1 = 0. 
Since TgTj is symmetric, it is sufficient to calculate the following components: 



N 



{T e Tj) 2k . 3 = cos jOi cosfcfl, 



A wiC-i) + w -i(»-i) w *d-i) + w -fc('-i) 

9 9 



AT 



= I^[ w o+fe)c-i) +w -o+fe)c-i) +w o- fe )('-!) +CL ,-o- fe )( ( - 1 )] 



4 

f (j = k) 

I o (j^k) • 

In a similar fashion, we have 

(Tot,; . . . : - o. 1 r„r,; L . y _, = | 



- J+ i - u, (i i e j 2 fc+ij+i - <; q (j ^ fc) 



Hence, we obtain (26). 

Next, we calculate the inverse of Tg. Let the matrix D satisfy 

T e T fl T = ^diag(2,l,...,l) = # 2 , (27) 



from which we can write 



D = y^diag(v / 2,1,...,1). 



Premultiplying and postmultiplying (27) by I? 1 respectively, we have 

D~ 1 T g Tg D^ 1 = I. 
We then know that D~ x Tq is an orthogonal matrix: 

D~ 1 Te(D~ 1 Te) T = I, 
from which we have (D~ 1 Tg)~ 1 = (D~ 1 Tg) T . Therefore, we obtain 

Tg 1 = TjD- 2 (28) 

with 

D- 2 := (D"Y = | diag(l !,...,!). 
We can obtain all the components of T ~ l explicitly from (28). □ 
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Remark 4. Arbitrary matrix Q can always be decomposed into 

Q = UP, 

where U denotes an orthogonal matrix and P a positive semidefinite symmetric 
matrix. This is said to be the polar decomposition of Q. 

We can redefine the sign of a part of the components of Tr and Tg as follows: 

(Tn)2k,2k = {T R )2k+\,2k+l = 7 ( 77- 

k \R 

{Tehkj = - cosfctf,-, {T e ) 2 k+i,j = - sinfcflj 

in Proposition 1. We see that K can be written as K = (TnD)(D~ 1 Tg) with 
the orthogonal matrix D~ 1 Tg and the diagonal matrix TrD whose diagonal 
components are positive. Hence, we know that K T = (D^ x Te) T (TrD) is the 
polar decomposition. 

In the following proposition, the determinant det(T0) with N = 2M + 1 is 
found in [2] and [7] . 

Proposition 3. Let N = 2M + 1. Then, the determinants of the matrices T# 
and Tg can be given by 

1 / R \ M ( M + 1 ) pjM+1/2 

det(Tfl)= (wUJ ' det(T ^^^- 

Hence, K = TrTq is not singular and K^ 1 can explicitly be written as 
{K- l h, = ^ (j = 1,2,..., N), 



(K -i, 2fc (R \ 



2k 


fRo 


~N 


KR 


2k 


(Ro 


iV 


\ R 



Proof. It is easy to obtain det(Tft) and det(Tg) by simple calculations: 



M 



det(Tfl) = 1 • H 



"k [r^ 



2 



1 f R\ 2Y -"= lk 



k=l 

1 / R \ M ( M+1 ) 



(Ml) 2 \R 



(Ml) 2 \RoJ 

and 



det(T e ) = ^dct^ 2 ) = JmrjrJ) = I ^ • (y 



iV\ 2M ] 1/2 N M+1 ' 2 



2 M 



From Proposition 2, we can immediately derive all the components of K 1 = 
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We notice that K is independent of the boundary Y and that K^ 1 is explicitly 
given, we can always directly obtain w from y without solving (24). 

The condition number of a matrix Q corresponding to the p-norm is defined 

as 

cond p (Q) = IIQIIpllQ- 1 ^, 

where || • || p is the matrix norm. When p = 2, the matrix norm HQH2 is the 
largest singular value of Q or the square root of the largest eigenvalue of Q T Q. 

Proposition 4. Let N = 2M + 1 and i? = R- Then, the condition numbers 
of T R , Tq and K = K 2 corresponding to the 2-norm are given as 

cond 2 (T R ) = M, cond 2 (T e ) = V2, cond 2 (K 2 ) = \/2M, 

which are independent of R. 

Proof. When we take Ro = R, we have 

T R = diag I 1,-1,-1,--,--,..., — -,-—]. 
& V ' ' 2' 2' M M J 

We can then easily see that ||Tb|| 2 = 1 and \\T^\\ 2 = M, which follows 
cond 2 (T R ) = M. 

We notice that the matrix norm HT0H2 is the square root of the largest 
eigenvalue of TgTj . Hence, from (26) in Proposition 2, we obtain \\Tg\\ 2 = y/N. 
Similarly, we have \\Tg \\ 2 = y/2/N. Therefore, we obtain cond 2 (T6i) = \/2. 

Next, since T R and TgTj are diagonal matrices, we can calculate 

K 2 Kl = {T R Tg){T R T e ) T = T R {T e Tj)T R 

_ N ( 1 1 1 I \ 

~ 2 H ' ' '^'^'■'■'^'MV ' 

Hence, we can see that H-K2II2 = VN and 1 1 J^T^ - 1 1 1 2 = V%M/\/N, which follows 
cond 2 (K) = \/2M. ' □ 

We assume that the boundary is a circle. Then, in a manner similar to the 
decomposition of K, the matrix S can be decomposed as follows: 

Proposition 5. Let p{9) be a constant p. Then, the matrix S can be decomposed 
into 

S = TgS Ro , 

where S Ro G jR (2M+i)x(2M+i) lg defined ag 

(S* )i,i = l, (S Ro )k,j = (k^j) 
{S Ra ) 2 k.2k = {S Ra ) 2 k+i.2k+\ = ( — ) {k = 1, 2, . . . , M). 
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We obtain the following proposition and theorem by similar discussions to 
Propositions 3 and 4: 

Proposition 6. Let N — 2M + 1. Under the assumption in Proposition 5, the 
determinants of the matrices Sr and S can be given by 



det(S Ra ) = 



Ro 
P 



M(M+1) 



det(S) 



N M+l/2 { M(M+1) 
2 M 



Hence, S is not singular. Namely, the approximate solution by the MTM is 
uniquely determined. 

Remark 5. We can solve (14) by Proposition 6. Using Propositions 2 and 5, 
y = S^(Tg')~ 1 f can be written in the form: 



a 



1 N 



N 



3=1 



ak = 



N \R 



k N 

cos kdj, 

3 = 1 
k N 



2 / \ 

bk= N\R ) E* sinfc ^ (* = 1,2,...,M), 

° 3=1 

which correspond to the approximation to the coefficients of the Fourier series 
(11): 



a 



i r 27r 



ak 



I ( Ji- 
lt \R 



k r 2ir 



f(9) cos k9 d6, 



6fc = K^) J^" W ahlk0d0 , (fe = l,2,...,M). 

Theorem 1. Under the assumption in Proposition 6, the condition number of 
S corresponding to the 2-norm is given as 



V2 



cond 2 (S') = < 



Ro 
P 

J_ /Ro 
I V2 V P 



_P_ 

Ro 

M-l 



M 

(Ro < p) 

( P <Ro< 2 1 /(2M) p) 

(2 l/(2M) p < Rq< x /2 p) 

M 

(Ro > V2p) 
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Proof. Since S Ro and TgTj arc diagonal matrices, we can calculate 



S T S 



Then, we know 
and 

\\s- 

Hence, we derive 




(TeS Ro ) - S Ro (T e Tj)S Ro 



Rq 
P 



\s\\ 2 , 



1 1 1 2 
1 1 2 



N 



max 



Rq 

P 



Rq 
P 



Ro 
p 



2M 



Rq 
P 



2M X 



- min J 2 I — ! j — 



Ro 
P 



2M 



2M 



-i -1 



\S\\o = 



and 



\s-% 



N 
2 

N_ 

2 

N 



N 
P 



P_ 

Ro 

JL 
Ro 

r 



M 



M 



Ro 
P 



(Ro < 2 1 /< 2M V) 
(Ro > 2 1 /( 2M )p) 

(Ro < P ) 

(p<Ro< V2p) 
(Ro > V2p) 



from which we obtain cond 2 (S') = ||5||2||<5 1 ||2- 



□ 



? opt 



We denote the characteristic length Ro that minimizes cond2(5*) by Rq 
which is called the optimal characteristic length. From Theorem 1, we obtain 
i?Q Pt as follows: 

Corollary 1. Under the assumption in Theorem 1, the optimal characteristic 
length is given by 

R o P t = 2 l/(2M) p _ 

Then, the minimal condition number of S is written as 

cond 2 (S) = 2< M - 1 )/( 2M ). 



(29) 



We can see that (29) is a monotonically increasing function of M and is less 
than Hence, the minimal condition number is less than cond2(>S) = V2 
for Rq = p. On the other hand, if we consider the conventional Trefftz method 
(Rq = 1) for the circular boundary whose radius is greater than 1 (p > 1), from 
Theorem 1 we have cond2(5) = \/2p M , which exponentially diverges for large 
M and p. Therefore, we know that the MTM is efficient. 

From Propositions 1 and 5, we immediately obtain the following theorem: 
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Theorem 2. Under the assumption in Theorem 1, the matrix SK corresponding 
to the MMFS can be decomposed into 

SK = TjKTg 

with the diagonal matrix A := Sr Tr € _r( 2M + 1 ) x ( 2M + 1 ) defined as 

A M = 1, Akj = (k f j) 
1 ( R\ k 

A 2 fc, 2 fe = A 2fc+1 , 2fc+ i = f -J (fc = 1, 2, . . . , M). 
Therefore, the approximate solution by the MMFS is uniquely determined. 

6. Numerical Experiments 

6.1. Condition number corresponding to the MTM 
We consider the following three domains (Figure 1): 

n k = {(r,e):r = pM(6), < 9 < 2^}, k= 1,2,3, 

where p^ (9) (0 < 6> < 27r) are defined as 

(a) circle: p { - 1 \9) = 1; {p m - m , Pmax) = (1,1) 

(b) ellipse: p^(9) = abj\J a 1 sin 2 6* + 6 2 cos 2 6> with a = 10,6= 5; (pmin, Pmax) = 
(5,10) 

(c) epitrochoid: p ( - 3) (9) = y/(a + b) 2 + 1 - 2(a + 6) cos(a9/b) with a = 3, 6 = 1; 

(PmimPmax) = (3,5) 

with p min = min p(6>) and p max = max ,c>(6>). 

0<e<27r O<0<2tt 




(a) (b) (c) 

Figure 1: Three domains 

We can confirm that the numerical results for k — 1 which are shown be- 
low are proper by comparing them with the corresponding theoretical results 
obtained in the previous section. 
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Figure 2 shows the condition number conc^S*) against < Ro < 15 for each 
domain tt k (k = 1, 2, 3) with N = 21 and M = 10. 




Figure 2: Condition number of S against Ro for each domain fij. (TV = 21, M = 10) 



From this result, we can observe that 
min cond2(5) = 

0<iJ <15 

We can see that the minimal cond2(5) and the corresponding Rq for k = 1 
coincide with the result of Corollary 1. When we use the conventional Trcfftz 
method (Rq = 1), the condition number cond2(S') for each domain ttk (k — 
1, 2, 3) with N = 21 and M — 10 can be obtain as follows: 

{1.4142 (jfe = 1) 

1.9430 x 10 9 (jfe = 2) . 
3.6112 x 10 6 (jfe = 3) 

Hence, we can see that the choice of a suitable Ro effects the drastic reduction 
of cond 2 (S'). 



1.366 at R = 1.035 (k = 1) 
167.439 at R = 5.850 (k = 2) 
39.481 at Ro = 3.636 (k = 3) 
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Figure 3: Optimal Rq against N for each domain Qf. 



Figure 3 shows i?Q Pt against N for each domain J7fc. In this figure, i?Q Pt 
seems to remain stable as N increases. We can observe that 

-Ro Pt ^ /°min = min p(6) 1 

which is theoretically true when k — 1 since R^ pt = 2 1 /( 2M ) j o = 2 1 /( JV ~ 1 )p is 
given by Corollary 1, and i?Q Pt ~ p — 1 for sufficiently large N. From this 
result, we can conclude that a suitable Rq should be a little bit greater than 

Pmin • 
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Figure 4: conc^S) for ijS pt against JV for each domain Qj. 

Figure 4 shows cond2(<S') for i?Q Pt against N for each domain fi^. Figure 5 
shows cond2(5') against N for each domain when Rq = 1. Comparing 
Figure 4 with Figure 5, we can confirm that the suitable characteristic length 
Rq can drastically reduce the condition number. 
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cond (K x ) (R=0 .5) 
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Figure 7: cond2(ii"i) (R = 0.5, 1.26) v.s. cond2(-fsT2) against N 

The matrix K depends only on (N, M, R, Rq). We notice that K is indepen- 
dent of fl. Figure 6 shows both of cond2(-?Ti) and cond2(i^2) against R when 
N = 21. From this figure, we can observe that cond2(ifi) has the minimal value 
2.22 at R = 1.26, which is smaller than cond 2 (.K2) = 14.1421 at any R. But, we 
cannot always take R that minimizes cond 2 (-ftTi) since the source radius R must 
be smaller than p m i n . Figure 7 shows both of cond 2 (ATi) for R — 0.5, 1.26 and 
cond2(if 2 ) against odd TV. We see from Figures 6 and 7 that cond 2 (if 2 ) is much 
smaller than cond 2 (-/fi) for most of i?, from which we confirm the efficiency of 
K 2 . 



6.3. Comparison of the MTM, the MFS and the MMFS 
We assume that the exact solution is given by 



u(x, y) = exp 



y 



+ y J \ % + y 

in the exterior domain Q — ^3 outside the epitrochoid boundary defined in the 
beginning of the section. We confirm the asymptotic behavior of u at infinity: 



u{x,y) 



exp 
exp 

1 H 



X' 4 
COS ( 

r 

cos 6 



y 



cos 



sm f 



= l + 0( 



r 



1 f sin 
1 - - 
2 



r 
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which does not satisfy the condition (3). Hence, we use the MFS to solve the 
exterior problem with respect to u := u — 1, which satisfies (3). Then, we can 
obtain the approximate solution to u = u + 1. 

We apply the following five methods for solving the problem: 

1. The MTM 

Solve Sy — f to substitute y into (11) 

2. The conventional MFS with the conventional basis functions (CMFS-CBF) 
Solve Aw = f to substitute w into (5) 

3. The conventional MFS with the modified basis functions (CMFS-MBF) 
Solve Aw — f to substitute w into (10) 

4. The modified MFS with the conventional basis functions (MMFS-CBF) 
Solve SK2W = f to substitute w into (5) 

5. The modified MFS with the modified basis functions (MMFS-MBF) 
Solve SK2W = f to substitute w into (10) 

We take (N,M) = (19,9). First, we confirm the accuracy of the solution 
obtained by all the methods. The absolute error between the numerical and the 
exact solutions on the circle whose radius is r is defined as 

e(r, 9) = \u(r cos 9, r sin#) — u(r cos 9, r sin#)|, 

where u stands for the numerical solution by one of the five methods. Figure 8 
shows the absolute errors e(10, 6) for R = 0.5, 1. We can see that the MMFS- 
MBF is the best accurate among all the methods. 




Figure 8: Absolute error {(R, N) = (1, 19), r = 0.5 (left), r = 1 (right)) 



Next, we confirm the accuracy of the solution obtained by all the methods 
in the whole exterior domain. We define the following maximum error: 

e(r) := max e(r, 9). 

0<6<2n 
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Figure 9: Maximum errors e(r) against the distance r from the origin 

Figure 9 shows the errors e(r) against r <E [10, 10 10 ] with (a) (R,N) = 
(0.002,7), (b) (R,N) = (0.2,11), (c) (R,N) = (0.5,19), (d) (R,N) = (1.2,19), 
(c) (R,N) = (1.2,11), (f) (R,N) = (1.0,19). We can observe from the figure 
that the both of the errors for the CMFS-CBF and the MMFS-CBF increase as 
r increases, while the both of the errors for the CMFS-MBF and the MMFS- 
MBF decrease as r increases, Hence, we know that the modified basis functions 
are suitable for approximating the solution in the whole exterior domain. 

Moreover, we can see from the same figure that in both of the cases where 
the conventional and the modified basis functions are used, the accuracy of the 
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MMFS is improved more than the accuracy of the conventional MFS. 

For small r, the MMFS and the modified basis function do not always give 
better accuracy. But, they give much better accuracy by taking suitable pa- 
rameters R and N. The accuracy of the MMFS-MBF is the best for the whole 
r in almost cases. 

Therefore, we can conclude that the MMFS-MBF is the best method among 
the five methods if we need to obtain a highly accurate solution in the whole 
domain. 

6.4- Condition numbers corresponding to the MFS and the MMFS 




Figure 10: cond2(A) v.s. cond2(S-ft") against R (N = 9) 
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Figure 11: cond2(j4) v.s. cond2(SK) against N (R = 1.2) 

Figure 10 compares the condition numbers of A and SK corresponding to 
the CMFS and the MMFS, respectively, when N — 9. Figure 11 compares 
the condition numbers of A and SK for R = 1.2, from which we can observe 
cond 2 (A) = 8.552 x 2.752^ and cond 2 (SJf) = 1.248 x 2.058^. Therefore, 
cond2(S 'K) is smaller than cond2(A). It is concluded that we can successfully 
reduce the condition number by using the MMFS. 

7. Conclusions 

We have proposed the MMFS with the modified basis functions for solving 
the exterior boundary value problem, based on the MTM. Under the assump- 
tion of the circular boundary, the condition number corresponding to the MTM 
is mathematically shown. Then, the optimal characteristic length Rq that mini- 
mizes the condition number is given in the mathematical form. The uniqueness 
of the approximate solutions by the MTM and the MMFS is also proven. 

The numerical experiments shows that the MMFS-MBF proposed in the 
paper is a more accurate method than the MTM, the CMFS-CBF and the 
CMFS-MBF. The condition number corresponding to the MMFS is smaller 
than the one corresponding to the CMFS. It is concluded that the MMFS-MBF 
is an efficient method. 

As future works, we will give some mathematical expressions of condition 
numbers that are not given in this paper. More methodical approach than one 
in this paper will be appeared in another paper. It is important to prove the 
reason why the MMFS can improve the accuracy. 
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